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Ion-exchange is commonly used to chemically strengthen glasses, by replacing small atoms by 
larger ones at sub T g temperature, thereby inducing a compressive stress. However, the resulting 
expansion of the glass remains lower than that predicted by the difference of molar volumes of the 
as-cooled glasses, an anomaly that remains poorly understood. Here, based on molecular dynamcis 
simulations of permanently densified sodium silicate glasses coupled with topological constraint 
theory, we show that the rigidity of the network controls the extent of the dilatation. Isostatic 
networks, which are rigid but free of eigenstress, show maximal expansion and, therefore, appear to 
be an attractive option to improve the toughness of glass. 

PACS numbers: 61.43.-j, 62.40.+i, 62.25.Mn 


PROBLEM DEFINITION: THE NEED FOR 
ULTRA-STRONG GLASSES 

Due to their non-crystallinity, glasses show unique 
properties in terms of optical, electronic and mechanical 
properties [I]. In particular, the transparency of usual 
silicate glasses is at the root of life changing applications 
like optical fibers or smartphones. However, even after 
decades of intensive research, glasses still break rather 
easily (smartphone drops often remain fatal) L2H3]. Their 
brittleness seriously limit the range of applications, which 
has been identified as a grand challenge for glasses M. 

Over the last decade, driven by the need for stronger 
screens for smartphones and tablets, the chemical 
strengthening of glasses through ion-exchange mm has 
emerged as the easiest and cheapest solution to this prob¬ 
lem. The most famous example is probably Corning® 
Gorilla® Glass |9l [10 . This process involves the replace¬ 
ment of smaller ions (typically Na + ) by larger ones (typ¬ 
ically K + ) at the surface of the glass, which is usually 
achieved by placing the glass in a salt bath (typically 
KNO 3 ) at a temperature below glass transition [IT]. This 
induces a compressive stress over tens of micrometers of 
the surface, which limits the risk of cracks propagating 
from surface flaws when subjected to tension. If poten¬ 
tial stress relaxation is neglected, the compressive stress 
profile obtained by ion-exchange follows: 


CT=-—(C(Z)- C „ g ) ( 1 ) 

when B is given by: 


B= l -I d P 

3 VdC 


( 2 ) 


where E is the Young’s modulus, v is the Poisson’s ra¬ 
tio, C(z) is the local concentration of substituted cation 


at a depth z from the surface, C avg is the average of 
such concentration over the whole range z- axis, and B is 
the linear network dilation coefficient, also known as the 
Cooper coefficient, in honor of Prof. Alfred R. Cooper 
m • The linear network dilation coefficient was defined 
in analogy to the thermal expansion coefficient with re¬ 
spect to the relative change of the molar volume V of the 
glass. Despite a large amount of empirical knowledge, 
the fact that B is typically lower than the value F? max 
that would be expected from the density difference be¬ 
tween the as-cooled alkali silicate glasses (e.g., between 
sodium and potassium silicate) remains poorly under¬ 
stood nn This is a serious limitation, as higher values of 
B would permit higher surface compressive stresses. This 
project specifically aims to better understand the role of 
the rigidity of the atomic network of silicate glasses on the 
relative amount of strengthening due to ion-exchange. 


METHOD: MOLECULAR DYNAMICS COUPLED 
WITH TOPOLOGICAL CONSTRAINTS 
ANALYSIS 


As accessing an atomistic picture of the effect of ion- 
exchange by conventional experimental techniques would 
be challenging, we rely on classical molecular dynamics 
(MD) simulations. Such simulations are analyzed within 
the framework of topological constraints theory (TCT) 
to characterize the rigidity of the network. TCT consid¬ 
ers the complex atomic networks of glasses analogous to 
simple mechanical trusses while filtering out the chemical 
details that ultimately do not affect macroscopic proper¬ 
ties Bang. The atoms undergo some mechanical con¬ 
straints arising from chemical interactions, which can be 
divided into two types: the 2 -body radial bond-stretching 
constraints and the 3-body angular bond-bending con¬ 
straints. Following Maxwell’s criterion of stability for 
mechanical trusses [16], glassy networks can either be 



2 




Pressure during cooling 


FIG. 2. Number of constraints per atom n c , for permanently 
densified sodium silicate glasses, with respect to the pressure 
applied during the cooling phase. The grey area indicates the 
domain of maximum expansion due to ion-exchange. 


RESULTS: COMPUTED ION-EXCHANGE 
COMPRESSIVE STRESS 


FIG. 1. Snapshot of the atomic structure of a sodium silicate 
glass cooled under zero pressure. Si, O, and Na atoms are 
represented in yellow, red, and blue, respectively. 


flexible, stressed-rigid, or isostatic, when the number of 
constraints per atom, n c , is lower, higher, or equal, re¬ 
spectively, to 3, the number of degrees of freedom per 
atom. 

To elucidate the influence of atomic rigidity on 
ion-exchange, sodium silicate glasses (Na 20 ) 3 o(Si 02 ) 7 o 
(noted NS hereafter, 3000 atoms, see Fig. [lj cooled un¬ 
der constant pressure P are prepared by MD, followed 
by relaxation to zero pressure at room temperature, but 
remaining permanently densified [IMS- As shown in 
Fig-i these glasses have been reported to show a rigid¬ 
ity transition with respect to P: flexible at low P and 
stressed-rigid at high P due to the change of the co¬ 
ordination number of Si and O atoms. This effectively 
delimits an isostatic pressure window around 8 GPa, in¬ 
side which NS shows a reversible glass transition [2QH24]. 
Typically, glasses show rigidity transitions according to 
composition, but a change of composition (e.g., the frac¬ 
tion of Na, or a replacement of Si by A1 atoms) would 
affect alkali atoms and therefore ion-exchange strength¬ 
ening, which would not allow us to identify the effect of 
network topology only. 


Taking the permanently densified NS systems as host 
glasses, we simulate the effect of ion-exchange by manu¬ 
ally replacing a given fraction of Na + by K + , following 
a methodology already well-established [T2l [25] . After 
ion-exchange, the glass is relaxed at zero pressure in or¬ 
der to compute the resulting increase of volume. All the 
processes are performed at 300K in order to avoid stress 
relaxation [26]. 

Fig. |3] shows the computed molar volume of ion- 
exchanged sodium silicate (cooled under zero pressure) 
with respect to the fraction of Na + replaced by K + ions. 
As suggested by Eq. |TJ the molar volume of the ion- 
exchanged glasses increases fairly linearly with the frac¬ 
tion of substituted cations. We note that, similar to 
what is observed in alkali aluminosilicate glasses [25] , 
once all Na + have been replaced by K + ions, the mo¬ 
lar volume of the ion-exchanged glass remains lower than 
that of the as-cooled potassium silicate glass (noted KS 
hereafter). This has been attributed to the fact that 
the ion-exchanged system is a forbidden glass m in 
which K atoms have a different coordination number 
than that in KS glass [12]. However, we observe that 
the molar volume of ion-exchanged NS glasses reaches 
97% of that of as-cooled KS glass, whereas, using the 
same methodology, the molar volume of ion-exchanged 
(Na 2 0)2o(Al20 3 )2o(Si02)60 glass reaches only 94% of 
that of the corresponding as-cooled potassium alumi¬ 
nosilicate glass [25] , This highlights the fact that net¬ 
work formers atoms and, more generally, the rigidity of 
the network influence the extent of volume expansion in- 
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FIG. 3. Molar volume of a sodium silicate glass, cooled 
under zero pressure, after replacement of a given fraction of 
Na + by K + ions. The red (blue) dashed line indicates the 
molar volume of as-cooled sodium (potassium) silicate. 


duced by ion-exchange. 

Following this idea, we perform the same analysis on 
permanently densified NS glass with different degrees of 
rigidity (see Fig. [2| by replacing 100% of Na + by K + 
in each case, and monitoring the resulting increase of 
molar volume. Fig. [4] shows the Cooper coefficient T>, as 
calculated from Eq. [2j with respect to the number of con¬ 
straints per atom n c before ion-exchange in the host NS 
glasses. Interestingly, B shows a broad maximum within 
n c = 3-3.4, which corresponds to a window of pressure 
during cooling of 8-16 GPa. This pressure window fea¬ 
tures a maximum of diffusion in the super-cooled liquid 
state m and has been shown to be characterized by an 
isostatic network (20]. We note that, in the pressure win¬ 
dow, the Cooper coefficient reaches its theoretical maxi¬ 
mum value, as obtained by the difference of molar volume 
between the as-cooled NS and KS glasses, which remains 
fairly constant with different cooling pressures. 


DISCUSSION: COMPETITION BETWEEN 
RIGIDITY AND STRESS 

Question remains about the origin of this maximum in 
the isostatic phase. We propose the following atomistic 
picture (as illustrated in Fig. [5]). (1) Flexible networks 
feature internal modes of deformation, or ffoppy modes 
[28) . which allow for local deformations with low energy 
barriers. The stress imposed by the stuffing of larger 
atoms can therefore be partially absorbed by local reor¬ 
ganizations of the network. (2) On the contrary, stressed- 
rigid networks are completely locked and cannot easily 
deform. Hence, rather than showing a global expansion, 
the induced stress is distributed as eigenstress [215] within 
the network where some of the bonds are stretched and 


FIG. 4. Cooper (linear network dilation) coefficient B af¬ 
ter substitution of all Na + by K + ions, with respect to the 
number of constraints per atom n c . The grey area indicates 
the domain of maximum expansion due to ion-exchange. The 
blue dashed line indicates B max , as obtained from the molar 
volume of as-cooled sodium and potassium silicate. The black 
dashed line serves as a guide for the eye. 


some others are compressed, as shown in Fig. E (3) 
Eventually, isostatic networks are characterized by a per¬ 
colation of rigidity [28] so that the displacement of one 
atom should affect the entire network, but are free of 
eigenstress. Consequently, they feature the largest ex¬ 
pansion upon ion-exchange. All together, this suggests 
that a better tuning of the rigidity of the host glass would 
allow us to reach higher compressive stress and, therefore, 
to design stronger glasses. 

FUTURE WORK 

For further research, we will try to find the structural 
signature of isostatic glasses to understand why they 
strengthen in an optimal way. This structural signature 
could be used as a metrics to quickly identify optimal 
compositions, e.g., from diffraction patterns. Although 
pure sodium silicate glass offers a simple starting point 
for us to better understand the underlying physics of ion- 
exchange, such a study will be extended to more realis¬ 
tic compositions of glass, e.g., modified aluminosilicate 
glasses. However, the rigidity of such glasses, so far un¬ 
known, will have to be evaluated first. Also, we keep in 
mind that higher surface compressive stress could also 
be compensated by higher stress relaxation. Although 
such relaxations are challenging to study directly via MD, 
they can potentially be predicted by accelerated relax¬ 
ation techniques m This option will be evaluated in 
the future. 

The authors acknowledge financial support for this re¬ 
search provisioned by the University of California, Los 
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FIG. 5. Schematic of the origin of the maximum compressive stress in isostatic networks. 


provisioned by the Physics of AmoRphous and Inorganic 
Solids Laboratory (PARISlab). 
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